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We review recent results on QCD at finite temperature. In particular we will discuss the chiral phase transition 
in two flavour QCD and new results on the excitation spectrum in the plasma phase of QCD. 
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1. INTRODUCTION 

One of the most exciting predictions of QCD is 
the existence of a phase transition at some criti- 
cal temperature to a new phase of strongly inter- 
acting matter - the quark gluon plasma (QGP). 
Lattice investigations of finite temperature QCD 
(FTQCD), which have been in progress for nearly 
15 years, have contributed a lot to our understan- 
ding of the phase transition and the nature of 
the high temperature phase. The determination 
of the phase transition temperature, the order 
of the phase transition and temperature depen- 
dence of bulk thermodynamic quantities are cen- 
tral targets of lattice investigations at finite tem- 
perature. Nonetheless, the analysis of FTQCD 
still leaves us with many open questions. This 
is, for instance, the case for most aspects of the 
phase transition in the presence of light quarks, 
for the electric screening of the heavy quark po- 
tential in the QCD plasma phase and even more 
so for the magnetic sector of QCD at high tempe- 
ratures. An understanding of the static as well as 
dynamic screening mechanisms in the high tem- 
perature phase is necessary in order to identify 
the fundamental excitation spectrum and quasi- 
particle structure of the QGP. All these questions 
have been addressed in lattice investigations over 
the past few years and the progress made, has 
been described in many excellent review talks at 
previous lattice conferences [1]. Still there are no 
final answers to the complex questions posed by 
FTQCD, and in this review we will again only be 
able to report on partial answers given by the nu- 
merous lattice studies performed during the past 
year. 



In the case of QCD with dynamical fermions 
the quantitative results on bulk thermodynamic 
quantities are by far not as detailed as in the pure 
gauge sector where subtle aspects of the phase 
transition like the scaling behaviour of the latent 
heat or the surface tension can now be investi- 
gated [2,3]. The simulations with dynamical fer- 
mions on large lattices have mainly concentrated 
on an understanding of the nature of the phase 
transition and its fiavour dependence. While it 
seems to be well established that the phase tran- 
sition to the high temperature chiral symmetric 
phase of QCD is first order in the case of four 
light fiavours [4] the situation is still unclear in 
the physically more interesting case of two light 
quark fiavours. We will discuss the present status 
of these investigations in the next section. 

A central goal of lattice investigations of the 
QGP is to understand the fundamental degrees of 
freedom in this phase and their interactions. Are 
these just weakly interacting quarks and gluons 
or do more complicated quasi-particle excitati- 
ons or perhaps even bound states with hadro- 
nic quantum numbers exist? An answer to this 
question can come through an analysis of spec- 
tral functions in the QGP and will require the 
understanding of the analytic structure of vari- 
ous correlation functions. During the past few 
years this program has been pursued by exami- 
ning, at finite temperature, the structure of corre- 
lation functions commonly used at zero tempera- 
ture to study the properties of hadrons (masses, 
wave functions). The progress made in under- 
standing the the analytic structure of such hadro- 
nic correlation functions as well as recent results 
for a related gluonic operator, the spatial string 



tension, will be discussed in section 3. In section 
4 we give our conclusions. 

2. TWO FLAVOUR QCD 

The QCD phase transition in the presence of 
light quarks is expected to be controlled by the 
chiral symmetry properties of the QCD Lagran- 
gian. In the limit of nf massless quarks the La- 
grangian has a global ?7^(1) x SU{nf) x SU{nf) 
symmetry. At zero temperature the axial Ua{^) 
is broken explicitly on the quantum level due 
to the non-vanishing anomaly while the chiral 
SU(nf) X SU(nf) flavour symmetry is broken 
spontaneously, giving rise to (n? — 1) massless 
Goldstone particles. Universality arguments for a 
generic effective Lagrangian reflecting these sym- 
metries suggest that for rif > 3 the QCD phase 
transition is first order, while in the case of nf = 2 
the situation is less clear [5]. If the ?7^(1) sym- 
metry is effectively restored at Tg, the transition 
is expected to be first order, while a second or- 
der transition can occur when the ?7^(1) symme- 
try only becomes realized at temperatures larger 
than Tg. The consequences of this latter scena- 
rio have been examined in some detail. Universal 
properties at the critical point are expected to be 
described by the critical exponents of the three 
dimensional 0(4) symmetry group [5,6] and an 
order parameter for the restoration of the ?7^(1) 
symmetry above Tg has been constructed [7]. 

The analysis of the chiral phase transition in 
two fiavour QCD on the lattice faces additional 
problems due the fact that neither in the stagge- 
red nor in the Wilson fermion formulation is the 
correct continuum chiral symmetry realized for 
finite lattice cut-off. The staggered fermion for- 
mulation does seem to have the advantage that 
at least a U(l) x U(l) subgroup of the conti- 
nuum symmetry group is realized. This insures 
that even for a finite cut-off a symmetry restoring 
phase transition takes place. However, its proper- 
ties might be infiuenced by the U(l) x U(l) sym- 
metry as long as one does not get close enough to 
the continuum limit. This would lead to a second 
order phase transition with 0(2) critical indices, 
which has explicitly been shown to be the case 
in the strong coupling limit [8]. The fact, that at 



present even in the simulations on the large 8x16^ 
and 8 x 32^ lattices with rather small quark mas- 
ses no first order phase transition has been found 
for the two fiavour theory with staggered fermi- 
ons [10,11] may thus not be the final answer for 
the two fiavour theory in the continuum limit. A 
careful analysis of the critical behaviour is nee- 
ded. 

In the Wilson fermion formulation there is no 
continuous symmetry left and a priori one would 
thus not expect to find a finite temperature phase 
transition at finite lattice cut-off ^ . Nonetheless 
a complicated phase structure has been found in 
the Wilson formulation as well [9,12,13] and its 
relation to the finite temperature chiral transition 
has to be understood. 

2.1. Staggered Fermions 

The simulations of two fiavour QCD on an 
8 X 16^ lattice with quark mass ma = 0.004, which 
have been performed by the Columbia group du- 
ring the past year [11,18] gave some indications 
for metastabilities in the time histories of the or- 
der parameter and the plaquette expectation va- 
lue at /3 = 5.48. This, however, did not find addi- 
tional support in the new simulations on a larger 
lattice of size 8 x 32^ at the same coupling, which 
have been presented at this conference [11]. These 
new data suggest that thermodynamic observa- 
bles take on values which are between the values 
previously identified as upper and lower edge for 
possible discontinuity at a first order transition 
on the 8 X 16^ lattice. Moreover, the average va- 
lues calculated on the 8 x 32^ lattice are in agree- 
ment with the average values extracted from the 
runs on the smaller lattice, if one averages over 
the whole data sample. This seems to rule out a 
first order phase transition. It should, however, 
be noted that the currently existing runs of the 
Columbia group on the large 8 x 32^ lattice are 
still rather short and the above conclusions thus 
have to be considered as being preliminary. 



We here talk about a true phase transition with non- 
analytic behaviour of thermodynamic quantities in the li- 
mit of infinite spatial volume. It should be stressed that 
such a transition exists for staggered fermions in the li- 
mit of vanishing quark mass and should not be confused 
with the pseudo-critical behaviour found in thermodyna- 
mic quantities for non vanishing values of the quark mass. 
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Table: Pseudo-critical couplings on various size lattices and for various values of the bare quark mass 
[10,14-20]. The numbers for Nr = 8 are based on an analysis of existing data for the chiral order 
parameter (i/''/') ^^ discussed in Ref. [21]. The last row gives the results from the fits shown in Fig. 2 



At present there are no indications for a first 
order transition in simulations with staggered fer- 
mions, which so far have been performed on lat- 
tice of size Nr X Nj^ with temporal extent Nr = 
4, 6, 8 and various values of the quark masses 
ma. On these lattices pseudo critical couplings 
l3c(Nr,ma) have been determined, which charac- 
terize the point of largest variation of the chiral 
order parameter, (xx); and/or the Polyakov loop 
expectation value. In the table we have collected 
the results for all pseudo critical couplings deter- 
mined so far for different parameter sets. As far 
as their dependence on the lattice size could be 
checked no systematic dependence has been found 
[19,21]. The results for the chiral order parameter 
on lattices of size 8 x N^ and three values of the 
quark mass are shown in Fig. 1. For ma = 0.0125 
this illustrates that there is no visible difference 
between the results obtained on lattices with spa- 
tial extent N^ = 12 and 16. 

Can we understand the quark mass dependence 
of the pseudo critical couplings in terms of sca- 
ling properties in the vicinity of a second order 
phase transition? Is the quark mass dependence 
compatible with the behaviour expected for a sy- 
stem with critical indices in the universality class 
of the three dimensional 0(N) spin models? In 
the presence of a symmetry breaking term in the 
action, as is the case for QCD with a non va- 
nishing quark mass, finite size scaling techniques 
are of little help to answer these questions. For 
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Figure 1. The chiral order parameter on lattices 
of size 8 X N^ and quark masses ma = 0.0125 
[16,17], ma = 0.00625 [10] and ma = 0.004 [18]. 
The results shown are for N^ = 16, except for the 
squares, which correspond to N^ = 12. 



sufficiently large lattices the correlation length in 
the system will be limited by the non vanishing 
quark masses rather than the size of the system. 
In order to study the scaling properties at the chi- 
ral phase transition of QCD with dynamical light 
quarks one thus has to follow an approach, which 
differs from the finite size scaling tests used in 



the case of the deconfinement transition in pure 
SU(N) gauge theories. We have to study the sca- 
ling properties in terms of the symmetry breaking 
fields (quark masses) and determine the scaling 
behaviour in the limit of vanishing quark mass 
[6,8,21]. The scaling properties can be deduced 
from the scaling behaviour of the singular part of 
the free energy 
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f,(t,h) = h-'f,(hyH,hy'^h) 



(1) 



where h is an arbitrary scale parameter, t = (T — 
Tc)/T denotes the reduced temperature and h is 
the coupling to the symmetry breaking field. The 
critical exponents j/j and yh are related to the 
more familiar exponents /3 and 6 through 
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In the continuum limit of lattice QCD the coup- 
lings t and h can be expressed in terms of the 
couplings appearing in the lattice Lagrangian, 
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h = maNj = niq/T . (3) 



The second derivatives of the free energy give 
the singular part of susceptibilities. Making use 
of the scaling properties of fg one finds for the 
scaling behaviour of the location of the peak of 
the susceptibilities [21] 
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+ c(maNry/'^'' . (4) 



Eq. 4 can be used to test the scaling behaviour of 
the pseudo critical couplings given in the table, 
although it should be stressed that in most cases 
these couplings have not been extracted directly 
through a calculation of susceptibilities but rather 
through an examination of the slope of the chiral 
order parameter. The result of such a comparison 
is shown in Fig. 2. A fit to the data set for the 
lattices with temporal extent Nr = 4 yields for 
the relevant combinationof critical exponents [21] 
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Figure 2. The pseudo critical couplings for two 
fiavour QCD as function of h = maNj and tem- 
poral lattice sizes Nj = 4, 6 and 8. 



which should be compared with the correspon- 
ding values 0.599 and 0.546 for the 0(2) and 0(4) 
symmetric models, respectively. Clearly, the pre- 
sently available data do not allow one to distin- 
guish between the two cases. However, even more 
remarkable than the proximity to the 0(N) ex- 
ponents is the fact that the shift of the pseudo 
critical couplings does indeed seem to depend on 
the dimensionless coupling h = maNj alone rat- 
her than on Nj and ma separately as it would be 
the case in the strong coupling limit. 

The presently available data for two fiavour 
QCD in the staggered formulation are consistent 
with the behaviour expected in the case of a se- 
cond order phase transition with 0(2) or 0(4) cri- 
tical indices. Moreover, there is no indication for 
a substantial difference in the critical behaviour 
between lattices with temporal extent Nj = 8, 
and the smaller Nj = 4 lattices. 

2.2. Wilson Fermions 

The analysis of the chiral phase transition with 
Wilson fermions is not as straightforward as it 
is with staggered fermions. To start with, it is 
not obvious that QCD with Wilson fermions will 
have a finite temperature phase transition on lat- 



tices with finite temporal extent Nr ■ Tliere is no 
global symmetry, which could be restored at high 
temperature, and thus one would not expect to 
find a second order phase transition in the infinite 
spatial volume limit of the lattice model. This si- 
tuation is quite different from pure gauge theories 
or QCD with staggered fermions, where the exi- 
sting global symmetries can lead to second order 
phase transitions at non vanishing lattice cut-off. 
If there is a finite temperature phase transition in 
the Wilson formulation at finite cut-off at all, one 
would thus expect to find only first order phase 
transitions. 

The strategy in studying the QCD finite tem- 
perature phase transition with Wilson fermions 
usually consists of two steps, which are analogous 
to the procedure used in the staggered fermion 
formulation: 

• On a lattice with given temporal extent 
Nr one identifies the pseudo critical va- 
lues of the hopping parameter, Kpc(Nr , g^), 
at which a crossover from low temperature 
thermodynamics to that of the high tem- 
perature phase can be observed. This may 
be done by studying the behaviour of the 
Polyakov loop, the pion screening mass ex- 
tracted from spatial correlation functions, 
or even by inspecting the behaviour of the 
plaquette expectation value. 

• One tries to follow this pseudo critical line 
to the point where Kpc reaches Kg. The lat- 
ter being defined as the zero temperature 
critical line, Kc(g^), where the pion mass 
vanishes at non vanishing lattice spacing. 

While the first step is pretty obvious, it is less 
clear whether the second step can be realized. 
What is the relevance of Kg for the thermodyna- 
mics on lattices with finite temporal extent, #^? 
Does it correspond to a critical point, which corre- 
ponds to a phase transition point with vanishing 
quark masses? Does a coupling 6/gl exist, where 
for a given value of Nr we find Kpc = Kg? 

The reason for the above questions already be- 
comes evident in the non-interacting (6/g^ -^ oo) 
weak coupling limit of lattice QCD with Wilson 
fermions. In this case the partition function on a 



Nr X Nj^ lattice is given by 

Z(k) = Y[[{l-8K + 4:Ksm'^(qo/2) 

3 

+4«^sm2(g,/2))' 

3 

-F4K2sin2go + 4K2^sin2g^]^ , (6) 

where the spatial momenta g^ are even multiples 
of ir/Na and go is non-zero and an odd multiple 
of ir/Nr due to the anti-periodic boundary condi- 
tions. For this reason one finds that even in the 
limit Na -^ oo the zeroes of the partition function 
remain in the complex plane. The zeroes closest 
to the origin are given by 



Kpc(Nr) 



8 - sin2(7r/2#^) ± isin(7r/#^) 
8 - sm^{Tr/2Nr)y + sm^{Tr/Nr) 



(7) 



and lead, in the limit Nr ^ cxd, to the familiar 
result Kc = 1/8. The free Wilson fermion theory 
thus has a singular point (zero of the partition 
function on the real axis) with a diverging corre- 
lation length only in the limit Nr -^ oo. For finite 
Nr, however, Eq. 7 defines only a set of pseudo 
critical couplings. In view of this it seems to be 
unclear whether a Kc(g^), at which, for instance, 
the pion screening mass (extracted from spatial 
correlation functions) becomes zero, exists for fi- 
nite Nr and non-zero g^ . 

The first simulations in which the finite tempe- 
rature transition with Wilson fermions has been 
studied [22] showed that the line of pseudo criti- 
cal values Kpc(Nr, g^) shifts into the strong coup- 
ling regime even for moderate values of the pion 
mass. The possibility has been raised that it 
would extend all the way down to 6/g^ = wit- 
hout leading to small values of the pion mass at 
the crossover point. The line of pseudo critical 
K-values has been examined by now in great de- 
tail for the case #^ = 4 [12,13,23,24]. In Fig. 3 
we show the result of a systematic study of the 
crossover behaviour in the Polyakov loop on lat- 
tices of size 4x8^ which has been performed by 
the MILC collaboration [13]. It can clearly be 
seen that the crossover is rather fiat for large and 
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Figure 3. Polyakov loop expectation value on l94V'9 
hopping parameter k [13]. 
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of size 4x8^ versus 6/g^ for various values of the 



small values of 6/g^ while there is a rather rapid 
crossover around 6/g^ = 4.9. A similar behaviour 
has been found for other observables. 

The QCDPAX collaboration has analyzed the 
behaviour of the pion screening mass on the 
crossover line. They find that this quantity is 
consistent with zero below 6/g^ = (3.9 — 4.0) and 
that the value determined for the pseudo criti- 
cal hopping parameter Kpc(4:,g^) in this region is 
consistent with the value Kg determined on lar- 
ger lattices from the vanishing of the pion mass 
[23,24]. 

In Fig. 4 we show a collection of values of the 
pseudo critical hopping parameter Kpc on lattices 
with temporal extent Nr = 4 and 6 and com- 
pare them with data for Kg. A collection of pion 
screening masses calculated from spatial correla- 
tion functions at Kpc as well as zero temperature 
pion masses calculated at the same k value is 
shown in Fig. 5. We note that this does not seem 
to suggest a vanishing of the screening masses in 
the vicinity of 6/g^ = 4.0, although Kpc clearly 
comes close to Kc in this regime. In any case we 
do not seem to have any indication of a true phase 
transition in this regime. 

While the general structure of the phase dia- 
gram found for Wilson fermions on lattices with 
Nr = 4 seems to confirm our expectation that 
there will be no genuine finite temperature phase 
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Figure 4. Pseudo critical hopping parameters de- 
termined on Nr = 4 and 6 lattices. The upper 
curve and the two fancy crosses indicate values 
for Kc determined at zero temperature. 



transition on lattices with finite temporal extent, 
the situation is complicated by the results found 
in simulations for Nr = 6 [25,26]. In this case a 
strong first order phase transition has been obser- 
ved for gauge couplings smaller than 6/g^ ~ 5.0. 
In fact, the results for Nr = 6 are drastically dif- 
ferent from those of the Nr = 4 simulations in 
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Figure 5. Screening masses extracted from spatial 
correlation functions on lattices with temporal ex- 
tent Nr = 4 and 6 and zero temperature pion 
masses determined on lattices with large tempo- 
ral extent. In both cases calculations have been 
performed at (or extrapolated to) the pseudo cri- 
tical hopping parameter value Kpc(Nr, Q/g^)- 



Figure 6. Plaquette expectation values 



Q/g' 



on various 



size lattices for 



0.18. In the 



low temperature phase the plaquette expectation 
values on Nr = 6 lattices agree with the zero tem- 
perature results up to the first order transition at 
l3 = 5.02. The vertical line indicates the location 

of Kr- 



this region. Some results for plaquette expecta- 
tion values are shown in Fig. 6. The differences 
in plaquette expectation values for Nr = 4 and 
6 as well as the gap at the first order transition 
for Nr = 6 (at l3 = 5.02 for k = 0.18) is much 
larger than what we know from simulations in 
the pure gauge sector or with staggered fermi- 
ons. This might be due to the fact that in the 
Wilson fermion case the simulations are perfor- 
med rather deep in the strong coupling regime. 
However, there also is another possibility, which 
at present cannot be ruled out: The large gap 
in the plaquette expectation values reminds one 
very strongly about the results obtained even de- 
eper in the strong coupling regime in simulati- 
ons with large number of fiavours [9,24]. In the 
case of nf > 6 strong first order phase transitions 
have been observed for finite Nr , which actually 
turned out to be temperature (Nr) independent. 
At present we cannot rule out that the first order 
transitions observed in the nj = 2 case on Nr = 6 
lattices are also temperature independent. This 



will certainly soon be clarified through simulati- 
ons with larger values of Nr ■ 

In the present simulations on lattices with 
Nr = 4 and 6 the crossover behaviour as well 
as the first order transitions occur at K-values at 
which the zero temperature pion masses are still 
quite large, i.e. to^/T>4.0. Following the line of 
pseudo critical couplings into the strong coupling 
regime does not seem to improve the situation. In 
fact, when comparing the zero temperature pion 
masses calculated at Kpc(Nr, g^) for Nr = 4 and 
6 it seems that below 6/g^ ~ 5.0 the pion mass 
in units of the temperature increases with incre- 
asing Nr- This is also shown in Fig. 5. It thus 
seems that one has to concentrate on the coup- 
ling regime 6/g^ > 5.0 for the analysis of the finite 
temperature phase transition with Wilson fermi- 
ons, if one tries to reach a regime where the pion 
mass is small in units of the temperature. Many 
of the peculiar results obtained with Wilson fer- 
mions so far seem to be due to specific features of 
the Wilson action in the strong coupling regime. 



It is the hope that the situation will improve clo- 
ser to the continuum limit. This, however, will 
require simulations with quite large values of Nr ■ 

3. THE HIGH TEMPERATURE PHASE 

Despite many years of analytical and numerical 
studies of finite temperature QCD little quantita- 
tive is known about the relevant excitations cha- 
racterizing the structure of hot QCD in the vici- 
nity of the phase transition to the QGP as well 
as in the high temperature phase itself. There 
are various indications that even at high tempe- 
ratures the excitation spectrum might be com- 
plicated in particular due to non-perturbative ef- 
fects in the magnetic sector of QCD. In order to 
make progress in our understanding of such non- 
perturbative effects the analytic structure of cor- 
relation functions has been studied in much detail 
over the last years. 

At finite temperature the temporal extent of 
the lattice is limited to the interval < r < l/T 
and thus correlation functions cannot be studied 
at large Euclidean time separations. Hence the 
studies focused on the analysis of long distance 
properties of spatial correlation function. A com- 
parison of the numerical results with model calcu- 
lations and high temperature perturbation theory 
then allows one to extract information on the ex- 
citation spectrum. 

During the past few years the structure of spa- 
tial correlation functions with hadronic quantum 
numbers has been examined in detail [27]. These 
studies have now also been extended to the ana- 
lysis of hadron propagators in temporal direction 
[28,29] and glueball operators [31]. A related cor- 
relator, wich actually has been studied several 
years ago [32], is the spatial Wilson loop. This 
simulation confirmed the prediction of an area 
law behaviour for spatial Wilson loops above Tg 
[33]. However, the detailed properties of the spa- 
tial string tension have only been examined this 
year [34-36]. We will discuss these new results in 
the next two subsections. 

3.1. The Spatial String Tension 

At high temperature an effective, three dimen- 
sional action can be derived by systematically in- 



tegrating out contributions from non-static mo- 
des in the QCD Lagrangian. The effective theory 
is a three dimensional Gauge-Higgs model with 
an adjoint Higgs fields resulting from the static 
temporal component of the gauge potential. The 
effective action takes on the form 



S, 



eff 



2#c 



{S^ + Sh + S,) , gl = gHT)T (8) 



where Su, denotes the action of a three dimensio- 
nal SU(Nc) gauge theory, Su is the Gauge-Higgs 
part of the action and Sy is a potential term for 
the Higgs field [37], which introduces in leading 
order quadratic and quartic self couplings for the 
Higgs field. The overall coupling, g^, is related to 
the temperature and the running coupling of the 
four dimensional theory. 

The string tension in this effective theory will, 
of course, be non- vanishing. Its dependence on gs 
and the coupling constants in Sy , however, might 
be complicated. The analysis of the spatial string 
tension in (3-|-l)-dimensional QCD at finite tem- 
perature thus will yield crucial information about 
the structure of the effective theory deduced in 
the framework of dimensional reduction. In fact, 
the early simulations suggested that the spatial 
string tension would be temperature independent 
or only weakly dependent on temperature. This 
would suggest a rather complicated structure for 
the dimensionally reduced theory where the Higgs 
sector and non-static modes would substantially 
infiuence the spatial string tension. Only recently 
spatial Wilson loops have been studied in the pure 
SU(2) gauge theory at temperatures substanti- 
ally larger than the critical temperature. Here it 
turns out that ag rises rapidly with T [34,35]. So- 
mewhat surprisingly the detailed analysis on large 
lattices and close to the continuum limit shows 
that the spatial string tension at large tempera- 
tures is even numerically close to what has been 
found for the string tension in three dimensio- 
nal SU(2) gauge theory [39]. One finds for the 
(3-|-l)-dimensional SU(2) gauge theory and tem- 
peratures T > 2Tc [40], 

\/cr,(T) = (0M9±0M4:)g\T)T , T>2Te .(9) 

where g^(T) is the 2-loop form of the tempera- 
ture dependent running coupling constant. Here 



the additional free parameter in the /3-function is 
determined as At = 0.0076(13)Tc. 

The spatial string tension has been extracted 
from pseudo potentials, Vt(R), constructed from 
space-like Wilson loops of size R x S 



W(R,S) 



(10) 



These pseudo potentials have been calculated on 
lattices of size Nr x 32^ at a gauge coupling, 
l3 = 2.74, which is known to be the critical 
coupling for the deconfinement phase transition 
on a lattice of size 16 x 32^ [38]. The latti- 
ces with temporal extent Nr at this value of the 
gauge coupling thus correspond to temperatures 
T/Tc = 16/ Nr- From the results for the pseudo 
potentials shown in Fig. 7 it is obvious that the 
slope at large distances, k(T), rises rapidly with 
T. This slope is related to ag through 



V<^siT)/T = y<f)#. 



(11) 



In Fig. 7b we show the ratio T j ^fal. It is ob- 
vious from this that a^ rises slower than linear 
in T . Note that the scale used for T jT^ in this 
figure is logarithmic. The linear rise seen in this 
figure thus suggests a logarithmic dependence of 
T I \J(Js{T) on T/Tc- In fact, the temperature de- 
pendence is well described as being proportional 
to the inverse running coupling constant g~^(T). 
A fit with such an ansatz yields the result given 
in Eq.9. 

This can be compared with the string tension 
of the three dimensional SU(2) gauge theory [39], 



(12) 



^= (0.3340 ±0.0025)(7^ 



In the gauge part of the effective theory obtained 
from dimensional reduction the coupling constant 
is given by g^ = g^(T)T. This result thus sug- 
gests that the spatial string tension is domina- 
ted by the pure gauge part of the effective three- 
dimensional theory. In fact, it has explicitly been 
checked that ag is rather insensitive to the struc- 
ture of the Higgs sector of Seff [41]. 

In the regime Tg < T < 2Tc the behaviour of 
(Ts differs from the above scaling behaviour. A 
fiux tube model [36,42] suggests that in this re- 
gime the rise of ag with T is due to the squeezing 



of the vibrating surface by the shorter temporal 
direction of the four dimensional lattice. The spa- 
tial string tension is then expected to behave as 



T 
<^siT) = — (7 



T < Tr<2Tr 



(13) 



This seems to be roughly fulfilled by the data 
shown in Fig. 7b. 

The general structure of ag discussed above has 
also been observed in much simpler gauge models 
like the Z(2) gauge model in (2-|-l) dimensions 
[36,42]. In this case one can also distinguish two 
regimes in the high temperature phase, in which 
(Ts shows a different functional dependence. In 
the region between Tg and 2Tc the rise in ag can 
be accounted for by the squeezing of the chromo- 
electric fiux tube due to the shorter temporal 
direction and is still characterized by the para- 
meters of the four dimensional theory, while for 
T > 2Tc fiuctuations are suppressed in the tem- 
poral direction and the behaviour of the string 
tension is controlled by the three dimensional ef- 
fective theory. 

Although the high temperature behaviour of 
the spatial string tension seems to be quite na- 
turally related to the temperature dependence of 
the pure gauge part of the effective three dimen- 
sional theory given by Eq. 8, the result is not 
entirely satisfactory. The temperature dependent 
running coupling constant, g^(T), extracted from 
the fits to (Ts turns out to be nearly a factor of 
two larger than the running coupling constant 
extracted from the heavy quark potential [43]. 
This alone is not disturbing, as different obser- 
vables may well lead to different values for the 
Ay-parameter appearing in g^(T) as long as the 
O(g^) corrections are not fixed. The fact that 
both quantities seem to be described by the ef- 
fective action, Eq. 8, with different choices of 
g^(T) is, however, unsatisfactory. It has been 
argued that the spatial string tension is not defi- 
ned through a static operator and a comparison 
with the string tension of the effective theory thus 
requires the inclusion of additional contributions 
from non-static parts of the Wilson loop operator 
[37,43]. In this case, however, it would be entirely 
unclear why the proportionality factor found for 
the spatial string tension is that close to the va- 
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(b) 



T/T„ 



Figure 7. Temperature dependence of the pseudo-potentials Vt(R) minus the constant self-energy con- 
tribution Vq on lattices of size Nr x 32^ for /3 = 2.74 versus the spatial separation R measured in lattice 
units (a). A fit to the linear part of these potentials yields the spatial string tension. The inverse square 
root of the spatial string tension in units of the temperature is shown in (b). 



lue found for the three dimensional theory. This 
point clearly requires further investigations. 

3.2. Hadron Correlators 

During the past few years the temperature de- 
pendence of correlation functions with hadronic 
quantum numbers has been studied intensively. 
Here attention has focused on the behaviour of 
two-point functions at large spatial separations. 
The aim has been to extract from the exponen- 
tial fall-off of these correlators screening mas- 
ses, which should yield information about the 
existence or non-existence of bound states with 
the quantum numbers selected by the correlation 
function used. In the baryon and (pseudo-)vector 
meson channels these screening masses turned out 
to be compatible with the propagation of nearly 
free quarks, possibly with some residual interac- 
tion among them. One finds 



niHadron - "Ifree - rnvT 



(14) 



with n = 2 for mesons and n = 3 for baryons, 
i.e. the screening masses are close to multiples of 
the lowest Matsubara frequency, po^min = ""T for 
free fermions, which characterizes the exponential 
fall-off of a fermion correlation function at large 
spatial separations. 



The (pseudo-)scalar sector, however, shows a 
markedly different behaviour. Here the screening 
masses are only about half as large as expec- 
ted from the propagation of two independent 
quarks. This might indicate that bound states 
with the quantum numbers of the pion still exist 
in the QCD plasma phase. A two-particle bosonic 
bound state would have a vanishing lowest Mat- 
subara frequency and thus could have a smaller 
screening mass. Some results for screening masses 
in four fiavour QCD obtained from spatial hadron 
correlation functions is summarized in Fig. 8. Si- 
milar results have been obtained for two fiavour 
as well as quenched QCD. 

In order to clarify the consequences for the exci- 
tation spectrum of QCD it is necessary to further 
study the analytic structure of the hadronic two- 
point functions. Recently the behaviour of these 
correlators at temporal separation [28,29] as well 
as their sensitivity to changes in the boundary 
conditions for the fermion fields has been studied 
[30]. Quite a different behaviour is expected in 
these cases as the screening masses will not be 
infiuenced by the non-vanishing Matsubara fre- 
quency. However, new difficulties in the interpre- 
tation of the correlation functions will show up. 
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Figure 8. Hadronic screening length in four- 
flavour QCD determined from spatial correlation 
functions on lattices of size 8 x 16^. 



At non-zero temperature hadronic correlation 
functions can only be studied in a short Eucli- 
dean time intervall, < r < 1/2T. As a con- 
sequence, higher momentum excitations will con- 
tribute to correlation functions, in which source 
and sink are seperated in the temporal direction. 
This will strongly affect the shape of the corre- 
lation function at all accessible distances, which 
can be seen explicitly in the case of a meson cor- 
relation function, constructed out of a quark and 
an anti-quark operator. To be speciflc we will in 
the following discuss meson correlation functions 
on Euclidean lattices of size Nr x N^ constructed 
within the staggered fermion formulation of lat- 
tice QCD. In particular we will describe numerical 
results obtained from simulations of four-flavour 
QCD using dynamical quarks of mass ma = 0.01 
on lattices of size 8 x 16^ [29]. The meson two- 
point functions can be constructed using local me- 
son operators, H(x), as sources. These are de- 
flned in terms of quark and anti-quark sources, 
X{x) and x(«), as 



H(x) = (t)H(x)x(x)x(x) 



(15) 



where (fi^ is the phase factor, which depends on 
the quantum numbers of H. We use correla- 
tion functions with separations in the a;o-direction 



(temporal direction), which project onto a flxed 
momentum, p, of the meson, 

G„ixo,p) = ^ J2 e'''^{H{xo,x)H{0,0)).{16) 



(xi,X2,X3) 

The correlation function, Gh, will receive contri- 
butions from quark flelds carrying arbitrary mo- 
menta, subject to the constraint that their sum 
is equal to the flxed external meson momentum. 
In the case of non-interacting quarks and anti- 
quarks the correlation function can be calculated 
explicitly. In particular one flnds for vanishing 
external meson momentum (p = 0) on an inflnite 
lattice 



Gh(xo,0) 



mr-^ 



^ y^TOsh^(y(2T^o^ ,.„, 

dy -2 .(17) 

cosh y 



This correlation function, together with results on 
flnite lattices with staggered fermions is shown in 
Fig. 9 on a logarithmic scale. The deviation from a 
simple exponential decay, which one would expect 
for a free meson correlation function, is clearly 
visible. In particular at short distances the high 
momentum modes from the free motion of the 
constituent quarks give rise to a steep rise of the 
correlation function. 





lo' 


1 ' 1 ' 1 ' 1 ' 


K 






8 


- 




10* 


r\ 


16 

-—32 


- 


H 


lo' 


^~~--!!S^ 


=o 


-: 






O 


10" 

lo' 


^^^^''^^""°^^=--==— _ 


- 




^^.--'-. \ \^''^ - 








■■-v<-^H\ 








10 


I.I.I. 







0.20 0.30 



Figure 9. Meson correlation functions in the 
high temperature limit (lowest order perturbation 
theory) on lattices with inflnite spatial extent and 
various values of Nj . Also shown is the result in 
the continuum limit {Nj -^ cxd) given by Eq. 17. 
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Figure 10. Pseudo-scalar (dots) and vector meson (squares) correlation functions constructed from point 
(upper row) and wall source operators (lower row) at two values of the gauge coupling in the QGP phase. 
Fig. (a) and (c) are for &/ g'^ = 5.3 whereas Fig. (b) and (d) correspond to &/ g'^ = 6.5. The open symbols 
denote the leading order perturbative result. 



In order to become sensitive to low momen- 
tum excitations in temporal correlation functions, 
which would allow one to determine the lowest 
mass in this quantum number channel, one eit- 
her has to perform fits with several exponential 
functions [28] or one has to construct correla- 
tion functions which project dominantly onto low 
momentum excitations. The latter approach has 
been followed in studies of the hadron spectrum 
at zero temperature, where one uses wall source 
operators. 






e.e 



'(i){x)x{xo,x+e)x{xo, a;-|-e').(18) 



Here e and e' denote all even sites in the 3- 
dimensional hyper-plane of the lattice. The per- 
turbative analysis of correlation functions with 
wall source operators shows that they only receive 



contributions from terms with small internal fer- 
mion momenta. In Fig. 10 we show results for 
temporal correlation functions using point sour- 
ces as well as wall sources. One clearly sees that 
with increasing value of the gauge coupling (in- 
creasing temperature) the wall source correlation 
functions become fiatter while the point source 
correlation functions become steeper. The latter 
refiects the importance of high momentum exci- 
tations in the point-point correlation functions. 
These higher excitations completely dominate the 
fall-off of the correlation functions at moderate 
distances and give rise to large effective masses. 
The wall source correlators on the other hand 
become sensitive to the low momentum excita- 
tions. In fact, the operator used in the correla- 
tion functions is constructed such that in the case 
of non-interacting constituent quark fields it pro- 



13 



jects onto the state where the constituent quarks 
carry vanishing momentum. 

In Fig. 11 we show results for the temporal mas- 
ses extracted from the wall source correlation 
functions and compare them with twice the ef- 
fective quark mass. As can be seen the temporal 
mass in the vector channel drops rapidly above Tg 
and is compatible with twice the effective quark 
mass. In the pseudo-scalar channel the tempo- 
ral mass also remains small. In fact, for T>Tc 
it still is smaller than twice the effective quark 
mass, which may indicate the existence of attrac- 
tive interactions in this channel. 

The results discussed here suggest that there 
are no bound states with large masses in the high 
temperature phase. In particular the screening 
masses extracted from the temporal correlation 
functions in the pseudo-scalar channel are a fac- 
tor two smaller than those found from the spa- 
tial correlators. This suggests that also in the 
pseudo-scalar channel the non-vanishing Matsub- 
ara frequencies for fermions strongly influences 
the screening mass in the pseudo-scalar channel. 
This gets further support from the analysis of spa- 
tial correlators in which the quark fields are forced 
to obey periodic boundary conditions in the tem- 
poral direction. Above Tg the screening masses 
extracted from such correlators are significantly 
smaller than those extracted from the correspon- 
ding operator with periodic boundary conditions 
for the fermion fields [29,30]. Below Tg, however, 
both type of correlation functions yield identical 
values for the screening masses. This indicates 
that the fundamental excitations are well locali- 
zed bosons in the low temperature phase. In the 
high temperature phase the fermionic substruc- 
ture becomes visible. 

3.3. Pseudo Wave Functions 

Another set of spatial correlation functions, 
whose structure has been examined during the 
past year are the so called spatial wave functions 
[44,45]. Similar to the spatial hadron correlation 
functions, these correlators are not directly rela- 
ted to the wave functions of hadrons. However, 
the fact that they still show an exponential decay 
[44], which is not the case in leading order pertur- 
bation theory, was again taken as indication for 
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Figure 11. Temporal screening masses in units of 
the temperature versus the gauge coupling /3 = 
Q/g^ for vector (squares) and pseudo-scalar (dots) 
mesons. Also shown is twice the effective quark 
mass in Landau gauge (triangles). The results 
shown are for four-fiavour QCD on a lattice of 
size 8 X 16^. The critical coupling for the chiral 
phase transition on this size lattice is /3c = 5.15. 



some non-perturbative effects in the high tempe- 
rature phase. It has been suggested that the be- 
haviour of these pseudo wave functions can also 
be understood in terms of the structure of the ef- 
fective three dimensional theory for the high tem- 
perature phase [45]. The pseudo wave functions 
would then describe the bound states of heavy 
quarks with a mass given by Eq. 14 (with n = 1), 
which are bound by a confining potential with 
a temperature dependent string tension, which 
at high temperatures is given by Eq. 9. Such a 
scenario suggests that the pseudo wave functions 
drop exponentially like 



|t/^(i?)| ~ exp (-V<T.(T)7rTi?3/2 



(19) 



The now existing quantitative results for the spa- 
tial string tension at high temperatures as well as 
the available information on the spatial screening 
masses will allow one to check the relation bet- 
ween the pseudo wave functions, the hadronic 
screening length and the spatial string tension in 
detail [46]. 
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4. CONCLUSIONS 



5. ACKNOWLEDGEMENTS 



In this talk we did not pay much attention 
to studies of thermodynamic observables like the 
transition temperature, the Debye screening mass 
or the temperature dependence of the equation of 
state. Of course, these quantities are of great im- 
portance for the phenomenology of the QGP and 
of direct relevance for the discussion of possible 
signatures for QGP formation in heavy ion expe- 
riments. The current status of numerical studies 
of theses observables has been discussed in detail 
in various reviews [1,47]. Certainly in these cases 
as well, we cannot be satisfied with the accuracy 
of presently existing results from lattice calculati- 
ons and more detailed studies of these quantitites 
will clearly be performed in the future. 

In this talk we tried to concentrate on some 
of the important conceptual questions in finite 
temperature QCD, which have to be answered 
with the help of numerical simulations. We have 
discussed the status of our understanding of the 
phase transition in two fiavour lattice QCD and 
the temperature dependence of various correla- 
tion functions in the high temperature phase of 
QCD. In the first case it is of great importance 
to reach a coherent description of the QCD phase 
transition in both commonly used fermion for- 
mulations. We thus have to pursue detailed in- 
vestigations of the phase structure of both theo- 
ries, although the calculations, which currently 
are performed close to the strong coupling regime, 
may not yet be of direct relevance for continuum 
physics. Similarly in the second case we have to 
study in more detail the analytic structure of cor- 
relation functions, which are not directly related 
to physical observables, with the aim to arrive 
at a consistent description of the excitation spec- 
trum of the high temperature phase QCD. 

We have, so far, not obtained a definite anwer 
from lattice simulations for any of these problems. 
Most likely we will have to wait for computers in 
the Terafiop range in order to come closer to a 
satisfactory answer. However, even the partial 
answers we have obtained so far are encouraging 
and important for the preparation of further in- 
vestigations along this line. 
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